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Abstract 

We explore the space of spherically symmetric, static solutions in the de- 
coupling limit of a class of non-linear covariant extensions of Fierz-Pauli mas- 
sive gravity obtained recently in I arXiv: 1007.0443] In general, several such 
solutions with various asymptotic limits exist. We find their approximate 
short and long-distance behaviour and use numerical analysis to match them 
at the Vainshtein radius, r*. Our findings indicate, that for a broad range 
of parameters, the theory does possess the Vainshtein mechanism, screening 
the scalar contribution to the gravitational force within r^,. In addition, there 
exists a class of solutions in the literature, for which the 1 /r gravitational po- 
tential is completely screened within the Vainshtein scale. However, numerical 
analysis indicates, that for this type of solutions, the gravitational potential 
does not decay at spatial infinity. 



1 Introduction and Summary 



It is a fundamental question of field theory, whether it is possible to write down 
a consistent, interacting lagranigian for a massive spin-2 state on Minkowski back- 
ground. The first attempt to answer this question dates back to the original work 
by Fierz and Pauli yj, who were able to find a unique ghost and tachyon - free 
mass term for a free Minkowski graviton. Breaking the gauge invariance of General 
Relativity (GR), the model propagates five degrees of freedom of a four-dimensional 
spin-2 particle with a nonzero mass m, as required by Poincare invariance. 

The Fierz-Pauli (FP) model however can not describe the real world. The prob- 
lem is connected with the observation of van Dam, Veltman and Zakharov (vDVZ) 
[21 E], that even a tiny mass of the graviton which one would expect to be important 
only in deep infrared, leads to an 0{1) modification of gravity at all scales. In- 
deed, in the m ^ limit, in addition to the usual GR interactions mediated by the 
helicity-2 component of the graviton, FP lagrangian describes the helicity-0 mode, 
coupled with gravitational strength to the trace of the external energy-momentum 
source, 

Crn^o ^ -lidnf + ^rcT. (1) 

Such a non-decouphng of the scalar degree of freedom would readily be ruled out 
by Solar System tests, which do not measure significant deviations from GR at the 
scales at hand. 

As first noticed by Vainshtein [1] however (see also [5]), in a generic nonlinear 
completion of the FP model, the linear approximation breaks down at a large dis- 
tance r^, from an astrophysical source, such that the solar system tests should be 
compared to the predictions of the theory in a highly nonlinear regime. The ef- 
fect originates in the strongly-coupled dynamics of the helicity-0 mode in a generic 
nonlinear massive gravity, which screens the scalar contribution to the gravitational 
potential. 

While restoring agreement with General Relativity, nonlinear dynamics of the 
helicity-0 mode leads to yet more severe problems, connected with the appearance 
of a ghost on a generic background, as first found by Boulware and Deser (BD) 
in the Hamiltonian formalism [S]. There, the ghost manifests itself in the loss of 
the Hamiltonian constraint, leading to propagation of an additional, sixth ghost-like 
degree of freedom on top of the five of the linear Fierz-Pauli theory. 

An alternative method of detecting the BD ghost - first proposed in [7J, is via the 
Stiickelberg treatment of the theory. In this approach, one restores the broken gauge 
invariance of massive gravity by introducing four scalars, encoding the longitudinal 
degrees of freedom of a massive graviton at high energies. The high-energy regime 
is then captured by going to the "decoupling" limit, in which only the operators 
suppressed by the lowest scale are retained. Schematically, the decoupling limit 
lagrangian for the helicity-0 graviton on Minkowski background can be written in 
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the following form for a generic extension of the FP model O El |9], 

^m-.o 3 -lidnf + ^^{d\f + ^ttT. (2) 

Here, A5 = {Mpui^Y^^ is the scale which is held finite in the limit at hand. As 
noted above, the cubic self-interaction of the helicity-0 mode screens its contribution 
to the gravitational potential of a source of mass M below the Vainshtein radius 

= {M/Mpin^y^^. Precisely due to this interaction however, the equation of 
motion for the scalar graviton is higher than quadratic in time derivatives, leading 
to the ill-posedness of the Cauchy problem and propagation of an extra scalar degree 
of freedom, which is inevitably a (BD) ghost. While infinitely heavy on Minkowski 
background, the ghost becomes light enough to destabilize any reasonable spherically 
symmetric solution in the theory - e.g. that of a localized lump of matter O HI |9]. 

Recently, a class of nonlinear generalizations of FP massive gravity has been 
proposed, which bypass these issues in the decoupling limit. By appropritely tuning 
the graviton potential, they cancel the dangerous d'^n- dependent operators order-by 
order in nonlinearity [TOKTT]. Moreover, the infinite number of contributions in the 
potential have been resummed into a two-parameter class of models, and the ghost- 
free property of the full theory has been proven up to and including the quartic order 
in nonlinearity in the Hamiltonian formalism [12] . The cosmological solutions in the 
decoupling limit of this class of models have been studied in [13] . Interestingly, this 
type of extensions of FP massive gravity can be obtained in the recently proposed 
framework of theories with auxiliary extra dimensions [T3| [T5] (for further studies 
of models with auxiliary dimensions, see [MHZ])- 

In the present letter we will explore the space of spherically symmetric, static so- 
lutions in the Minkowski-space decoupling limit of the class of ghost-free extensions 
of the Fierz-Pauli model, found in [lOl [HI [12]. In this limit, one recovers a two- 
parameter scalar-tensor theory of a very special form, determined by symmetries at 
hand. We will call this class of models "A3 - theories" below. The decoupling limit 
action is interesting in several ways; it represents a unique theory, to which any 
nonlinear, ghost-free massive gravity should reduce at high energies. Moreover, it 
is a certain scalar-tensor generalization (see [11]) of higher- derivative, yet ghost-free 
scalar interactions, the simplest of which - the cubic one, was found in the decou- 
pling limit of the DGP model [18] in ref. [19]. The cubic helicity-0 self- interaction 
of the DGP graviton has been extended to incorporate the full set of terms with 
similar properties - the so-called Galileons ^U\. Galileons are well known to enjoy 
special properties, such as second-order equations of motion, nonrenormalization of 
lowest-order interactions, presence of the Vainshtein mechanism, etc. - all of these 
properties (including the Vainshtein mechanism, as we find below) are characteristic 
to the A3 theories as well. Finally, the decoupling limit theory admits a ghost-free 
self-accelerated solution, at the same time providing a novel mechanism of decou- 
pling the scalar degree of freedom from arbitrary matter [13]. 

The study of an analogous problem has recently been conducted in [22], [231 [21] • 
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It has been found by Koyama, Niz and Tasinato [23], that for a class of parameters of 
the theory (described by vanishing /3 - a particular parameter, to be defined below), 
the theory admits the Vainshtein mechanism, successfully restoring agreement with 
GR to a very high precision at sub- Vainshtein scales. Moreover, for /3 7^ 0, a solution 
within the Vainshtein radius has been found in [21], for which the 1/r contribution 
to the gravitational potential is completely screened, making it observationally un- 
acceptable. In the present note we further study the space of spherically symmetric, 
static solutions in this class of theories and show that even if /3 7^ 0, for a broad 
subset of parameters (defined by /3 < 0), another solution exists, which does exhibit 
the Vainshtein mechanism. Furthermore, we also find that the solution of [21] with 
the screened 1/r potential within the Vainshtein radius, matches a non-decaying 
solution outside. 

The letter is organized as follows. In section 2 we describe the framework and 
some technical details, involved in the discussion. We derive the basic equations and 
comment on their similarity with those in Galileon theories. Section 3 deals with 
the discussion of spherically symmetric, static solutions. We find that in general, 
the theory admits multiple such solutions in the small and large distance regimes. 
Finally, we investigate their matching using numerical analysis. 

2 The Formalism and Technical Analysis 

Our starting point is the two-parameter decoupling limit action of the ghost-free 
class of generalizations of Fierz-Pauli theory [TO] [TT] , 

^ = -lh>^'£::u'ha, + h^-'X^^J + ^^h^'Xj^J + ^^h^-'Xj^J + -^h^'T,.. (3) 

Here, the first term represents the usual GR kinetic term for the tensor mode, a and 
f3 are two arbitrary constants and A3 = {Mprn^y^^ is the scale that is held fixed 
in the limit at hand; by {£h)^^ we denote the linearized Einstein operator acting 
on metric perturbation, S'^^h^p = —\{Pih^u — d^dah^ — d^dah"^ + d^dyh — ri^i,Dh + 

The three identically conserved symmetric tensors X^"^[n] depend on second 
derivatives of the helicity-0 field 11^;^ = d^d^n in the following way, 

X^^^ = --£ "'"'e ^ U ^ 

^^fiu 2 pa'-'-alJ-t 

^ 1^ ap7p /3<^ TT TT 
liv 2 "I'-'-afi'-'-pa-i 

xS) = e/''%/'^'n„;,n,.n^5. 

Being irrelevant for the present analysis, we ignore the helicity-1 mode of the massive 
graviton everywhere below; it enters the full nonlinear action at least quadratically, 
therefore allowing for a consistent solution for which it is not excited at all. 
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Before turning to the discussion of spherically symmetric solutions, let us briefly 
comment on symmetries of the action Under linearized diffeomorphisms, the 
helicity-2 graviton transforms in the usual way, 

K. ^ K. + d,iu + d^i,, (4) 

while the scalar mode tt is invariant in the decoupling limitQ. The diff invariance 
is then manifest in the fact that the currents X*^") are identically conserved. On 
the other hand, there is also a manifest Galilean symmetry, under which the scalar 
mode shifts as 

TT — )■ TT + x^h^ + c, h^^ = invariant (5) 

with bn and c constants. Galilean transformations however can be redefined in a way, 
that realizes this symmetry only up to a total derivative, [TTj. Indeed, by making a 
helicity-2 redefinition, h^^ h'^^ + rj^^Tc + ad^irduiT / A.^ , one can eliminate the hX^^^ 

and hX^'^^ couplings in (jS]), while the hX^^^ coupling can not be eliminatecil; as a 
result, one generates the full set of Galileon self-interactions for the scalar mode, as 
well as an additional coupling of the form S^vrS^vrT^'^ on top of the usual coupling to 
the trace of the external energy-momentum source Then, even if we redefine the 
transformations so that h'^^ is invariant, the Galilean invariance will still persist 
up to a total derivative - precisely due to the Galileon structure of the decoupling 
limit theory in the Einstein frame. 

All of these facts indicate, that the spherically symmetric solutions in the given 
class of theories should not differ dramatically from those in modified gravities a la 
Galileon. In particular, for the subspace of parameters defined by /3 = 0, for which 
the theory reduces to Galileon-modified gravity, this statement is automatic and 
massive gravity therefore possesses all well-known nice properties of Galileon field 
theories - the presence of Vainshtein mechanism among them. For /3 7^ on the 
other hand, the situation is more subtle, and we will discuss this case in great detail 
below. 

Varying the action with respect to /i^" yields, 

while the equation of motion for the scalar mode reads as follows, 

8 8ahf^''( — -F '^P'^F ^ -\- —F "'P'^F '^'^ TT + S— F "'P'^F I^^^T\ T\ A — {] 

^This should not come as a surprise, since the helcity-0 field represents a physical longitudinal 
polarization of the massive graviton in the decoupling limit. 

^The "coupling" is of course the usual conformal mixing of the helicity-2 and -0 modes, 

its elimination corresponding to transition to Einstein frame. 

•^This coupling is irrelevant for a static, non-relativistic source. However it may play a significant 
role for gravitational lensing for instance, |25j. 
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The latter equation can be rewritten in terms of the hnearized curvature invariants 
in the following way, 

_^(iw)n^^n,;3 = o, (7) 

where R^^J , G^^J = E'^^hafi and R^^^^p denote respectively the linearized Ricci, Ein- 
stein and Riemann tensors composed of the helicity-2 field h^^. Presence of the 
Riemann tensor in the latter expression makes it impossible in general to alge- 
braically reduce the system of equations of motion to a set of pure - vr and pure - 
h^^ equations, whereas for /3 = this is possible. This is just another manifestation 

of the fact found in [IT] , that one can eliminate the h^'^XiiJ couplings by a nonlin- 

(3) 

ear redefinition of the tensor mode, while the coupling h^^X^iu , if present, can not 
be absorbed by any such redefinition of variables. 

The most general spherically symmetric ansatz for the metric perturbation is 
given as follows, 

hoo = a(r), hij = f{r)Sij + b{r)ninj, (8) 

where a, f and b are general functions of the radial variable, while rii denotes the 
unit vector in the radial direction. We will exploit the diff invariance of ([3]) to set 
b{r) = 0. Substituting the ansatz into the 0-0 component of the Einstein's equation 
([n]), one obtains, 

;^,(-rV' + rV-^.M=-|(.f)'^^r„. (9) 

where a prime denotes a derivative by r. The spatial components of (E]) on the other 
hand, reduce to the following two equations, 

(^rf — ra' 
(iL-sL 

\t r 

These equations are consistent only if the integration constants vanish, leading to a 
single equation, 

rf -ra' -2rn' + -^{Ti'f = 0. (10) 

^3 

Finally, the equation of motion for the helicity-0 field vr is of the form d^J^ = 
{r'^J^)'/r'^ = 0, where denotes the conserved Noether current of the shift sym- 
metry vr ^ vr + c. Since the equations of motion that follow from the lagrangian 
([3]) contain no more than two derivatives per field, can not contain higher than 
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a single derivative of vr, leading to an algebraic equation for vr', just as in the case 
of Galileons [211] • Plugging the anzats ([H]) into eq. ([7]), the scalar equation reads, 



2rV' - r'a' + 2^rvr' {a' - f) + Q^a'{^'f) = 0- (H) 

Below we explore the space of spherically symmetric, static solutions to the 
system (|9l)- (|TT]) . We will be interested in a spherical source of size constant 
density p and negligible pressure, so that the only nonvanishing component of the 
energy-momentum tensor is Too = P d{R). Integrating eqs. ([9]) and ffTTl) subject 
to boundary conditions a(0) = /(O) = 7r(0) = 0, one reduces the system to the 
following set of equations for the radial derivatives of the fields. 



a' = -^-A^rA(l + 2/3A^), (12) 



M(r) 



/' = -174 + ^3^^(1 - - (13) 



-W fi^x2,^9^2 omj^s To;o2x5 _ / (v)^(l + 6/3A^) Outside the source . . 



[1 + 6/3A ) Inside the source 



where M(r) denotes the mass of the source within a sphere of radius r, and we have 
defined, 

^ = TJZ^ ■ (15) 



3 Solutions 

We will first concentrate on solutions to the system (IT2l)- (fT^ outside the source. 
The helicity-0 equation f lT^ is quintic in A; it can not be solved exactly and in 
general we should expect multiple solutions. Outside the source, there exists a dis- 
tance scale - the Vainshtein radius r^,, at which a typical solution changes regime; we 
will therefore divide the space further into regions outside and inside the Vainshtein 
radius and solve the system f |T2l) -( |T^ in each of these regions separately. The non- 
trivial part is to find an appropriate matching of solutions in the two regions, which 
in principle requires an exact solution to the system. Having no handle on exact 
analytic solutions, we will implement numerical analysis in support of our findings. 
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3.1 Solutions beyond the Vainshtein radius 

For r ^ r*, there exist two types of solutions. One of them (which we call "Solution 
1" below) is obtained by simply neglecting all nonlinearities in A in eq (1141) . so that 

with corrections suppressed by powers oir^/r. This is an asymptotically flat solution 
of the linear FP theory - fields vanish at infinity, and it manifestly exhibits the vDVZ 
discontinuity. This type of asymtpotic behaviour one encounters in straightforward 
nonlinear generalizations of the FP masisve gravity, see e.g. [261 EZl EH [29] . 

Since the right hand side of (|T4l) contains a small overall factor (r*/r)^, formally 
a second solution (" Solution 2" ) is possible, for which the left hand side vanishes in 
the zeroth approximation, 

A oc vr'/r = C ^ TT ~ A^r^, a, / ~ A^r^, (17) 

with C a nonzero constant and up to corrections suppressed by powers of r*/r. This 
is an asymptotically non-decaying solution - both the helicity-0 and helicity-2 fields 
increase with distance on this branch (see also the discussion in section 3.3); the 
fact of its existence will play an important role below. 

3.2 Solutions within the Vainshtein radius 

Well within the Vainshtein radius, r*/r ^ 1, and for /3 < 0, there exists a solution 
("Solution 3") for which the r.h.s. of eq. (fT^ vanishes in the zeroth approximation. 
On this branch, the leading r-depencence of the fields reads as follows. 




TT ~ A^r^, a = f = + O 



M 




and one recovers with great precision the gravitational potential of GR, up to tiny 
corrections due to the helicity-0 graviton. One should note, that these corrections are 
smaller than analogous ones in the DGP model (see [301 El])- For /3 < therefore, 
the theory admits the Vainshtein mechanism! 

Finally, there exists an interesting solution found by Koyama, Tasinato and Niz 
("Solution 4") [22| |21]. for which the usual 1/r part of the gravitational potential is 
completely screened inside the Vainshtein radius. Mathematically, this corresponds 
to last terms on each side of eq. ( fT4l) being dominant over the rest of the contribu- 
tions. Then, the leading dependence on r in a and / cancels, leaving the following 
expressions for the solution at hand, 

( M M 

V 2/3; r ' Mprl ' ^ Mpr, ^ ^ 
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This modifies the gravitational potential within to an observationally unaccept- 
able form. 

Within the source, as it follows from (I14p . A is constant, matching the value at 
r = Rof a. particular solution inside the Vainshtein radius discussed above (Solution 

3 or 4). This means, that vr and its first derivative are continuous across the surface 
of the source. 

3.3 Matching at the Vainshtein radius 

Matching the solutions found inside and outside the Vainshtein radius is perhaps 
the most nontrivial task to perform. As remarked above, this procedure can not be 
carried out analytically and numerical analysis should be addressed. We therefore 
will not be able to claim full generality in this section - we will rather present results 
for a particular choice of the parameters a and /3; these results seem pretty generic 
- at least to the extent that numerical analysis can provide. Our findings indicate, 
that if /3 < 0, Vainshtein mechanism works (at least for all values of the parameters 
consistent with this condition, that we have been able to check). This means that 
the GR solution (Solution 3) which exhibits the Vainshtein mechanism inside r^,, 
does match the asymptotically flat Solution 1, exhibiting the vDVZ discontinuity 
beyond the Vainshtein scale. One example of this is shown in fig.l, which gives the 
dependence of A on r/r* for a = —1 and /3 = —1. The asymptotics of the plotted 
solution at r/r* ^0 and r/r* — )■ oo clearly coincide with the Solution 3 ( A ~ const) 
and Solution 1 (A oc respectively. On the other hand, the asymptotically non- 
decaying Solution 2 (A ~ const) outside the Vainshtein radius matches the Solution 

4 (A oc 1/r) inside, as shown on an analogous plot in fig. 2 for the same values of 
the parameters. For /3 > 0, only the latter type of solutions exists (see fig. 3). 

It is important to address the fact, that the leading terms in the expressions 
for A for Solns. 2 and 3 do not depend on the mass of the source. One can then 
ask the question: how physical are these solutions, since the fields do not vanish 
even in the M — )■ limit? In this limit, the non-decaying Solution 2 describes an 
asymptotically non-fiat configuration with nonzero (and increasing with distance) 
gravitational potential throughout the whole space. On the other hand, the region 
of validity of the the sub- Vainshtein Solution 3, r <^ r*, shrinks to zero for the 
vanishing mass of the source. Moreover, for any nonzero M, corrections to the 
gravitational potential are tiny within the Vainshtein radius, as noted above @ . 
Finally, for M = 0, there of course exists a physical solution, on which all fields are 
identically zero. 

^In fact, these arguments are valid only if one does not consider quantum corrections. In the 
full quantum theory, one can make predictions only down to distances of order of the inverse cutoff 
of the effective theory, Vq ^ ^3^- For Planckian sources, rq merges with the Vainshtein radius r* 
and the theory transitions from the linear regime straight into the quantum one, in which definite 
predictions can not be made. For any transplanckian source M 3> Mp on the other hand, the 
Solution 3 describes a well defined region in which the helicity-0 contribution is screened and one 
recovers GR with great precision, yet quantum corrections are under control. 
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A comment on yet another class of solutions on asymptotically non-flat back- 
grounds may be in place here; on a self-accelerated background of , the hX^^^ 
mixing vanishes (while the kinetic term for the helicity-0 graviton comes from the 
curvature of the backround spacetime), making it possible for the scalar mode vr not 
to be excited at all by a localized source, thus decoupling it from arbitrary matter 
(for details, see [13j). 
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Figure 1: The solution exhibiting the Vainshtein mechanism for a 
/3 = -l. 
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Figure 2: The asymptotically non-decaying solution for a = —1 and /3 = — 1. 



In summary, our findings indicate that in general there are multiple spherically 
symmetric, static solutions in the Minkowski-space decoupling limit of ghost-free 
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Figure 3: The asymptotically non-decaying solution for a = —1 and (3 = 1. No 
asymptotically decaying solutions exist for this choice of parameters. 



nonlinear extensions of the Fierz-Pauli model, found in [TT | [T2]. Moreover, for a 
wide range of parameters, defined by /3 < 0, the Vainshtein mechanism seems to work 
in the theory - there exists an asymptotically-flat solution, which screens the con- 
tribution of the helicity-0 mode at sub- Vainshtein scales, successfully hiding it from 
the Solar System tests. In addition, there exists an asymptotically non-decaying 
solution, which totally screens the 1/r potential inside the Vainshtein radius |24] . 
dramatically reducing the gravitational attraction at these scales - this scenario is 
clearly ruled out by observations. 
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